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INTRODUCTION 

In current and advanced turbofan engines, increased cycle pressure and tempera- 
ture are employed to achieve reduced specific fuel consumption and increased 
thrust/weight ratios. As a result, turbine airfoils are subjected to increasingly 
higher heat loads which escalate the cooling requirements in order to satisfy life 
goals. If high turbine efficiency is to be achieved, cooling flow requirements must 
be kept as low as as possible and the anticipated significant effect on the local 
metal temperatures (and ultimately on airfoil life) must be evaluated. One way to 
keep the quantity of cooling air bounded would be to develop more efficient internal 
cooling schemes. Also, development of reliable methods for predicting the thermal 
and aerodynamic performance of these schemes would eliminate the necessity of main- 
taining a safety margin surplus of cooling air. 

To predict local metal temperatures, and ultimately blade life, the turbine de- 
signer faces the problem of accurately evaluating the local heat transfer from the 
hot gas to the blades and from the blades to the coolant flowing within the internal 
passages of the airfoil. Furthermore, for an acceptable design, the heat transfer 
must be tailored to reduce gradients while at the same time maintaining an appro- 
priate overall level of metal temperature. To accomplish these design goals, the de- 
signer must be able to predict local boundary conditions both on the external hot 
gas side and the internal coolant side of the airfoil walls. These boundary condi- 
tions, which consist of the external gas driving temperature (Tg), external heat 
transfer coefficient (h ex ), internal heat transfer coefficient (he), and internal 
coolant temperature (Tc), are input to computer prediction codes that solve the heat 
conduction-structural aspect of the airfoil analysis, subsequently arriving at a 
life prediction. 

To predict local coolant side heat transfer coefficients and coolant temperature 
rise in the cooling passages of current and advanced turbine blades is difficult be- 
cause of the complex geometry of the passages and the turbulence-promoting devices 
that are used to achieve high levels of airfoil cooling effectiveness. Currently the 
design analysis of airfoil internal passage heat transfer and pressure drop relies 
on correlations derived for the most part from testing models in a static 
(non -rotating) environment. Executing tests with rotation is difficult and costly 
and, as a consequence, there are almost no data in the literature that a turbine de- 
signer can use with confidence to account for the effects of rotation on the in- 
ternal heat transfer coefficients and pressure loss in typical turbine blade de- 
signs. Some data are available for smooth tubes over a limited range of relevant pa- 
rameters, but application of these data to the complicated flow passages of a tur- 
bine airfoil would not be appropriate. As a consequence, adjustment factors are gen- 
erally applied to the static test derived correlations to bring them into nominal 
correspondence with engine experience. This, in theory, accounts for rotation ef- 
fects. Design application of computer codes to predict 3D viscous passage flow has 
been limited because the computer codes are in an early stage of evolution. These 
codes, however, offer the potential for analyzing arbitrary geometries accounting 
for all real-world effects. 
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OBJECTIVE 


The objective of this 36-month experimental and analytical program is to develop 
a heat transfer and pressure drop database, computational fluid dynamic techniques 
and correlations for multi-pass rotating coolant passages with and without flow tur- 
bulators. The experimental effort will be focused on the simulation of configura- 
tions and conditions expected in the blades of advanced aircraft high pressure tur- 
bines so that the effects of Coriolis and buoyancy forces on the coolant side flow 
can be rationally included in the design of turbine blades. 

EXPERIMENTAL MODEL 

The coolant passage heat transfer model features a four-pass serpentine arrange- 
ment designed to reflect the passages in a gas turbine blade. A photograph of the 
partially assembled model is shown in figure 1. The radial passages have a square 
cross section with 12.7 by 12.7 mm (0.5 by 0.5 in.) dimensions. The model test sec- 
tion surfaces consist of copper segments, 10.7 by 49.3 mm (0.42 by 1.94 in.), which 
are heated with thin film electrical heaters. The copper segments are separated from 
each other with 1.5 mm (0.060 in.) strips of G-ll composite material. The copper 
segments are held in place with a frame constructed of G-10 composite material. G-10 
and G-ll were chosen because of the low thermal conductivity and high strength. 

The heat transfer coefficients on the test section surfaces are obtained from 
(a) heat balances on the copper segments and (b) thin film heat flux gages mounted 
on the segments. Pressure measurements are obtained along the straight sections and 
in the turns. The model instrumentation with 64 heated segments and 16 pressure 
measurement locations, is shown in figure 2. The 12 heat flux gages are positioned 
to determine local heat transfer gradients in the inlet region and on the walls 
downstream of the first two turns.' The electrical leads from the heaters, thermo- 
couples and heat flux gages are brought through hermetically-sealed connectors and 
sliprings to a computer-based data acquisition system. The pneumatic leads are 
brought through a leak-tight connector to 16 differential pressure transducers (ZOC 
model by Scani valve) and a reference pressure transducer. The electrical signals 
from these pressure transducers are also routed through the sliprings to the data 
acquisition system. 


DATA REDUCTION 

Data acquisition /analysis consists of three general categories: equipment cali- 
bration, model heat loss measurement, and heat transfer coefficient calculations. 
The equipment calibration follows standard experimental procedures. 

Next, model heat losses are determined experimentally. The model is brought up 
to a constant wall temperature, steady-state operating condition with no coolant 
flow. The energy needed to maintain constant wall temperatures then equals the 
energy conducted through the model walls to the environment. These values are 
normalized by the temperature difference between the heated copper walls and the 
model exterior to account for slight variations in environmental conditions during 
testing. This difference represents the driving gradient for conduction losses and 
is measured for each test. 
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For each copper element the net energy convected to the fluid is calculated by 
subtracting the electrical line losses and conducted heat losses from the total 
energy supplied. Bulk fluid temperatures are then calculated based on an energy bal- 
ance for each flow path section as follows: 


T 


b 


out 


^net, 4 walls 
itc 

P 



where the model inlet bulk temperature is measured. Once bulk fluid temperatures are 
determined, heat transfer coefficients are calculated from the equation 


q net, wall 


where T[> is the average of the inlet and exit bulk temperatures. Thus heat trans- 
fer coefficients are calculated for each individual copper wall element. This quan- 
tifies local heat transfer processes and helps isolate local effects of rotation. 


RESULTS 


To date experimental heat transfer and pressure drop results have been obtained 
for stationary conditions (no model rotation). This data will serve as the baseline 
for assessing the impact of rotation and buoyancy effects. Figure 3 presents the 
heat transfer results obtained for a Reynolds number of 25,000 which is a represen- 
tative nominal value for airfoil cooling passages. This figure shows that the heat 
transfer downstream of the entrance and turns decays monotonically to approximately 
the level (marked by arrows) for fully developed flow in a square duct as estab- 
lished by Lowdermilk, et al. (ref. 1). The expected periodicity from pass-to-pass is 
achieved. 


Heat transfer in the turns is approximately 2 to 2 1/2 times that of fully de- 
veloped flow levels. There is no heat transfer data in the literature with which to 
directly compare the current results. Sahm and Metzger (ref. 2) have found heat 
transfer levels around 2 to 2 1/2 times fully developed values for 180-degree square 
cornered bends. This is in agreement with the current results. Figure 4 compares 
the entry length behavior for the current data with results from the literature. For 
the first outward straight passage the four data points (solid circles) include the 
guard heater region and the following three straight segments. As the flow enters 
the channel from the plenum it passes through a screen and begins boundary layer de- 
velopment. Previous measurements just downstream of the screen showed the velocity 
profile of the core flow to be similar to fully developed turbulent flow. Because 
the near wall shear flow characteristics would be similar to those of developing 
flow, it would be expected that the entry length heat transfer behavior be in rea- 
sonable agreement with results for combined hydrodynamic and thermal boundary layer 
development. This is supported by the comparisons of figure 4 (Boelter et al, 
ref. 3, and Aladyev, ref. 4 for pipe flow; Yang and Liao, ref. 5, for square duct 
flow). Heat transfer levels downstream of the first bend would be expected to be 
lower than for the first passage. This, again, is in agreement with observation. A 
creditable correlation for the effects of the bend secondary flows on the downstream 
heat transfer is not available. 
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Pressure drop measurements are shown in figure 5 for a non-heated stationary 
test condition at the same Reynolds number as for the heat transfer results but at a 
lower pressure. The results indicate a negligible pressure drop in the straight 
sections which is in agreement with predictions and about 0.8 dynamic heads pressure 
drop through each turn. This agrees with a value of 0.71 dynamic heads determined 
from reference 6 for connected 90-degree square elbows. The slightly higher loss 
for the model is believed to result from the diffusions and accelerations the flow 
experiences as it negotiates the area variations in the turns. 

THREE-DIMENSIONAL VISCOUS FLOW COMPUTATIONS 

Computation of heat transfer and pressure drop for the reported experimental 
conditions employing a generalized aerothermal fluid dynamic solver for three- 
dimensional elliptic, turbulent, steady flows has been started; results are 
currently not available. 

TESTING WITH ROTATION 

Table I defines the planned test points for the first phase (smooth passage 
investigation) of this program and reflects the domain of interest for turbine blade 
coolant passages. 
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TABLE I 


Test Matrix 



— 

Bas ic 



Secondary Dimension- 



Dimensionless 

Parameters 

less Parameters 


Test 

R 


AT_ 

H 

Ap ftH 

Gr 

G r max 

Comments 

Mo. 

e 

lx o 

Tin 

d 

P v 

Re^ 

X 10“ 8 

1 

25,000 

0 

0.14 

— 

0 

0 

0 


2 

12,500 

0 

0.14 

— 

0 

0 

0 

r Nonrotating 

3 

50,000 

0 

0.14 

— 

0 

0 

0 

) 

4 

25,000 

0.210 

0.14 

52 

1.53 

0.32 

2.01 

Baseline 

5 

12,500 

0.210 

0.14 

52 

1.53 

0.32 

0.50 

l R e Varied 

6 

50,000 

0.210 

0.14 

52 

1.53 

0.32 

8.03 

S 

7 

25,000 

0.420 

0.14 

52 

3.06 

1.28 

8.03 

[ R 0 Varied 

8 

25,000 

0.105 

0.14 

52 

0.76 

0.08 

0.50 

i 

9 

25,000 

0.210 

0.07 

52 

1.53 

0.16 

1.00 

l AT/T Varied 

10 

25,000 

0.210 

0.21 

52 

1.53 

0.48 

3.01 

( 

■ 11 1 

25,000 

0.210 

0.14 

34 

1.00 

0.21 

1.31 

H/d Varied 

12 

25,000 

0.210 

0.14 

52 

1.53 

0.32 

2.01 

l a Varied 

13 

25,000 

0.210 

0.14 

52 

1.53 

0.32 

2.01 



Dimensionless Parameters 


Reynolds number, Re 

Rotation number, Ro 

Density or temperature 
ratio 

Radius ratio 

Rotational Grashof 
number , Gr 


pVd H /p = Vd H /v = md H /(Ay) 
Sld H /V 

Ap/p, AT/T 
H/d H 

Ro 2 Re 2 (AT/T)(H/d H ) 
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PHOTOGRAPH OF UNINSTRUMENTED COOLANT PASSAGE HEAT TRANSFER MODEL 


TRAILING EDGE ( + «) PLANE TEST SECTIONS REMOVED 


ORIGINAL PAGE IS 
OF POOR QUALITY 



Figure 1 


84 - 7 - 37-2 


ORIGINAL PAGE IS 
OF POOR QUALITY 

INSTRUMENTATION PLAN FOR COOLANT PASSAGE HEAT TRANSFER MODEL 


TEST SECTION ELEMENT IDENTIFICATION: 

SURFACES 1-32 ARE ON SIDE WALLS PERPENDICULAR TO VIEW SHOWN 
SURFACES 33-48 ARE ON +8 LEADING PLANE 
SURFACES (49M64) ARE ON " + ST TRAILING PLANE 
PRESSURE MEASUREMENT LOCATIONS |~1~| - ~ 

NOTE: EACH TEST SECTION SURFACE IS INSTRUMENTED WITH TWO THERMOCOUPLES 
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NUSSELT NUMBER, Nu f = hd/k. 


0 = 0.0 rpm 
P = 146,9 psia 

m = 0,01 31 Ib/sec 


HEAT TRANSFER RESULTS 


= 25260 


T wall = 618 ° R 

T b- inlet = 538 ° R 
AT /T = 0.149 


Re f, inlet 

R o = °-° 
H/d h = _ 


LEADING AND TRAILING TEST SECTION SURFACES 



SIDE WALL TEST SECTION SURFACES 



Fi gure 3 
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WALL SECTION NUSSELT NUMBER 


2.50 


o 


Re f, inlet = 25,260 
Nu f = 0.01 74 Ref 0 8 


YANG & LIAO: SQUARE DUCT 
MEASURED DATA Re - 23.500 


ORIGINAL PAGE IS 
OF POOR QUALITY 


BOELTER, YOUNG & IVERSEN: 
COMBINED HYDRODYNAMIC 
+ THERMAL DEVELOPMENT 
Re - 50.000 
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ALADYEV: CIRCULAR DUCT EXPERIMENT 
WITH WATER, Re - 20,000 


# FIRST OUTWARD STRAIGHT 
SEGMENT - 4 WALL AVERAGE 



Figure 4 


Entry Length 


Heat Transfer Comparison 


PRESSURE PROP RESULTS 


U o rprn 
P = 36 2 psia 
m = 0 0123 Ib/sec 


•lT = 0 0 # F 
Re = 25045 


V in ^ 39 3 It/sec 

Q'P = 4 38 psl = 0 843 .n H ? G 



Figure 5 


409 


